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O ■■ Abstract. 



We summarize recent significant progress in the development of a first-principles 
Qj ' formalism to describe the formation and evolution of matter in very high energy heavy 

i-Q . ion collisions. The key role of quantum fluctuations both before and after a collision is 

emphasized. Systematic computations are now feasible to address early time dynamics 
CN ' essential for quantifying properties of strongly interacting quark-gluon matter. 

> 

CN ■ PACS numbers 

O ! 1. Introduction 

The evolution of matter produced in the "little bang" of heavy ion collisions, to borrow 
a term from a recent cosmological model [1] , is a process of Ekpyrosi4j|, in which highly 
occupied coherent gluon states in the incoming nuclei are excited, producing very hot 
and dense matter shortly after the collision. The analog of the decaying inflaton field 
in this conflagration is a highly occupied Glasma [21 [Sj H] field. Quantum fluctuations 
play an essential role: a dominant class of these are responsible for the energy evolution 
of the Color Glass Condensate (CGC) nuclear states before the collision. The other 
dominant fluctuations grow exponentially with the square root of the proper time after 
the collision and likely isotropize matter in the Glasma. The mechanism of isotropization 
and possible thermalization bears close analogy to the inflationary dynamics of "pre- 
heating" in the early universe [5l[6l|7]. We will outline here recent work quantifying the 
role of quantum fluctuations in the early time dynamics in heavy ion collisions 

I A term from Stoic cosmogony wherein the universe is destroyed in a fiery conflagration and 
subsequently reconstituted at regular intervals. Heavy ion collisions containing billions of events 
repeatedly undergo this process and are therefore closer in spirit to the Ekpyrotic scenario than the 
origin of the universe in an initial singularity. 
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2. Factorization and resummation of quantum fluctuations 

The treatment of high energy nuclear colhsions in quantum field theory is motivated by 
the natural separation of scales between large x static color sources pi 2 in the nuclear 
wavefunctions and small x dynamical gauge fields A^'". The evolution of this separation 
of scales leads to a Wilsonian renormalization group description formalized as the CGC 
effective field theory [8]. At high energies, the color sources are strong {pi 2 ~ Vfi') where 
(7 ^ 1 is the QCD coupling) and color fields become time dependent immediately after 
the collision. Inclusive quantities, computed in the presence of these time dependent 
color fieldd§l, are expressed in terms of retarded propagators, thereby allowing real time 
computations with initial data at negative infinity [9| [TO]. 

Before the collision, the effective ground states of the nuclei are coherent states 
described as classical fields with occupation numbers 0(1/(7^). In QCD, a classical 
description of Ekpyrosis is obtained by solving Yang-Mills equations with static light 
front color sources [21 [HI [12]. At this classical level, energy dependence is introduced 
"by hand" via the saturation scale Qs, which is closely related to the variance of the 
Gaussian random distribution of color charges in the MV model [HI [H] of high nuclear 
wavefunctions. The energy dependence arises at next-to-leading (NLO) order from 
quantum fiuctuations (of relative strength g) about the classical background of the 
nuclei. These naively sub-leading contributions are enhanced by logarithms, which give 
a5ln(xi2) ~ 1 for small 0:1^2 and have to be resummed to all orders in perturbation 
theory. In the strong field regime, there is an additional resummation {gpi^2)"' at 
each order in the as expansion. These respectively radiative and multiple scattering 
contributions are generated by the JIMWLK renormalization group (RG) equation, 
QiJ^^ ) = "^1,2 Vr[pi,2], where Vr[pi,2] are gauge invariant weight functionals for the 
distributions of sources in the two nuclei, and ^{1^2 is the corresponding JIMWLK 
Hamiltonian [151 [IB] • 

An important consideration is to show that quantum fiuctuations from the two 
nuclei don't talk to each other before the collision. Albeit required by causality, it is 
by no means assured in a weak coupling treatment; we are able to show formally that 
factorization of the contributions of the weight functionals Vr[pi^2] to inclusive quantities 
is obtained at leading logs in x accuracy [HI [TBI [IE]. These quantum modes are boost- 
invariant p^ = modes, where p"^ is Fourier conjugate to the space-time rapidity. 

In the process of Ekpyrosis, nuclear coherence is lost and p^ 7^ modes are 
generated. These modes are generically unstable [201 1211 [221 [22], and grow in an 
expanding system as (asexp(2y//Ir))", with p ~ Qs, where n denotes the order in 
perturbation theory beyond the classical leading order contribution. These "leading 
instabilities" are comparable to the background field at r ~ I/Q5 and have to be 
resummed to all orders, leading to qualitatively different behavior. 

This early time behavior has a robust analogy to "pre-heating" dynamics in 

§ Examples of similar dynamics include Schwinger's mechanism for electron-positron pair production 
in strong QED fields and Hawking radiation in the vicinity of the event horizon of a black hole. 
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inflation |S1 El [7] where quantum fluctuations, enhanced by parametric resonance with 
the background Inflaton fleld, are conjectured to play an important role in "turbulent 
thermalization" of the early universe. In particular, quantum evolution of such a system 
can be represented by the stochastic average over repeated evolution of the inflaton fleld 
to which quantum seeds drawn from a Gaussian distributed ensemble of fluctuations 
are added [211 [25]. As we shall discuss, this formulations bears an exact analogy to the 
treatment of the problem in QCD [26] . 

3. Results and Outlook 

After factorization and resummation of leading logs in xi^2 and leading instabilities 
asexp(2y/JIf), the energy-momentum tensor is expressed as 

{TnLL.+Linst.= J[DpiDp2] W^^W^M J[Va] Fo[a] TI^^[A[pi, p2]+a]{x) .(1) 

The weight functionals H^xi,2[Pi,2] satisfy the JIMWLK RG equation discussed 
previously. The argument A = {A, E) denotes collectively the components of the 
classical flelds and their canonically conjugate momenta on the initial proper time 
surface; analytical expressions for these are available at r = 0''' [21 [3]. The initial 
spectrum of fluctuations Fo[a], Gaussian in the quantum fluctuations a, has a variance 
given by the small fluctuation propagator in the Glasma background fleld as r — )■ 0+. 
In practice, the path integral in a is determined by solving the classical Yang-Mills 
equations repeatedly with the initial conditions at r = 0+ given by 

Ainit. = Anit. + J dp^ [c^ a^ (x) + cl a^* (x)] . (2) 

Here A collectively denotes the quantum flelds and their canonical conjugate momenta. 
The coefficients c^^, with K collectively denoting the quantum numbers labeling the 
basis of solutions, are random Gaussian-distributed complex numbers given by 

^kC) = ^5..., , k'^.'> = k-C) = 0. (3) 



9 KK' ' \ J^ K' / \ K K' 

Explicit expressions for the small fluctuations and their conjugate momenta, denoted 
here by a^{x) were obtained in ref. [26]. The inner product of these solutions satisfles the 
orthogonality condition (a^|a^^, ) = N^ 5^^, with the measure dp^^ (a mix of integrals 
and discrete sums) that ensures J dp^ N^^ 6^^, = 1. 

In ref. [26], a numerical algorithm was outlined to compute eq. ([I]), thereby 
describing both Ekpyrosis and inflationary dynamics including essential leading 
quantum fluctuations. (We emphasize the formalism holds for any inclusive quantity in 
heavy-ion collisions including for instance parton energy loss and sphaleron transitions 
at early times.) We can now study how lumpy initial conditions for color charges 
Pi 2 in nuclear wavefunctions at the energy /rapidity of interest transform into the flow 
of matter. The early universe analogy is becoming more robust experimentally with a 
"WMAP-like" spectrum of spatial anisotropics now being resolved in heavy ion data |27j . 

We carried out an extensive study of the formalism outlined here for a scalar 0^ 
model which, among several QCD-like features, has a spectrum of unstable quantum 
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modes which are amphfied by resonant interactions with the background field [22]. In (j)'^ 
(and other scale invariant theories), the amplitude of the field is inversely proportional 
to its period of oscillation. Slightly different amplitudes, corresponding to different 
quantum seeds, lead to differing oscillation periods; a stochastic average over these 
leads to decoherence in the evolution. A striking consequence is hydrodynamic flow 
with an ideal equation of state. Eq. (|2]) is a realization of Berry's conjecture [29] which 
is believed to be necessary for thermalization [30] of a quantum system. Thermalization 
and onset of quasi-particle dynamics have been studied in the scalar theory; similar 
studies are feasible in the QCD framework of ref. [26] . 
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